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We review some basic notions and results of White Noise Analysis that are used in 
the construction of the Feynman integrand as a generalized White Noise functional. We 
show that the Feynman integrand for the harmonic oscillator in an external potential is 
' a Hida distribution. 
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1 Introduction 
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Path integrals are a useful tool in many branches of theoretical physics including quantum me- 
chanics, quantum field theory and polymer physics. We are interested in a rigorous treatment 
of such path integrals. As our basic example we think of a quantum mechanical particle. 

On one hand it is possible to represent solutions of the heat equation by a path integral 
representation, based on the Wiener measure in a mathematically rigorous way. This is 
stated by the famous Feynman Kac formula. On the other hand there have been a lot of 
attempts to write solutions of the Schrodinger equation as a Feynman (path) integral in a 
useful mathematical sense. The methods are always more involved and less direct than in 
the euclidean (i.e. Feynman Kac) case. Among them are analytic continuation, limits of 
finite dimensional approximations and Fourier transform. We are not interested in giving full 
reference on various theories of Feynman integrals (a brief survey can be found in [2]) but we 
like to mention the method in [1] using Fresnel integrals. Here we have chosen a white noise 
approach. 

White noise analysis is a mathematical framework which offers various generalizations 
of concepts known from finite dimensional analysis, among them are differential operators 
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and Fourier transform. Although we will give a brief introduction to white noise calculus in 
section 2 the reader unfamiliar with this topic is recommended to the monographs [5], [15], 
[4] and the introductory articles [12], [17], [19], [21]. 

The idea of realizing Feynman integrals within the White Noise framework goes back to 
[6] . The " average over all paths" is performed with a Hida distribution as the weight (instead 
of a measure). The existence of such Hida distributions corresponding to Feynman integrands 
has been established in [3]. In [8] Khandekar and Streit moved beyond the existence theorem 
by giving an explicit construction for a large class of potentials including singular ones. 
Basically they constructed a strong Dyson series for the Feynman integrand in the space of 
Hida Distributions. This approach only works for one space dimension. Their construction 
was generalized to (one dimensional) time-dependent potentials of noncompact support in 
[13]. 

In this work, which originated in the White Noise workshop on Madeira in August 1993, 
we carry those ideas over to perturbations of the harmonic oscillator. Hence instead of 
constructing a Dyson series around the free particle Feynman integrand we expand around 
the Feynman integrand of the harmonic oscillator as obtained in [3] . The external potentials 
to which the oscillator is submitted correspond to the wide class of time-dependent singular 
potentials treated in [13]. 

In [1, chap 5] the path integral of the anharmonic oscillator is defined within the theory 
of Fresnel integrals. Compared to our ansatz this procedure has the advantage of being 
manifestly independent of the space dimension. Despite the lack of a generalization to higher 
dimensional quantum systems our construction has some interesting features: 

• The admissible potentials may be very singular. 

• We are not restricted to smooth initial wave functions and may thus study the propa- 
gator directly. 

• Instead of giving a meaning to the Feynman integral we define the Feynman integrand as 
a Hida distribution. By taking expectation we get the propagator. On the other hand 
one may now use the toolbox of white noise analysis and apply differential operators to 
derive variational relations or Ehrenfest's theorem, see [5, chap 12], [18]. 

2 White Noise Analysis 

The starting-point of White Noise Analysis is the real Gel'fand triple 

S (R) C L 2 (R) C S' (R) , 

where S' (R) denotes the real Schwartz space. Using Minlos' theorem we construct the White 
Noise measure space (S' (R) ,B,/j) by fixing the characteristic functional in the following way: 

C(f) = J expi(u;,f) ^H = exp(-i^/ 2 (r) dr) , / € S (R) . 

We denote by (-, ■) the bilinear pairing between S' (R) and S (R) and by |-| the norm on 
L 2 (R). 
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Within this formalism a version of Wiener's Brownian motion is given by 

B (t) := (u, l [0 ,t)) = J^uj{s) ds . 

We now consider the space (L 2 ), which is defined to be the complex Hilbert space 
I? {S' (R) , B, fi). For applications the space (L 2 ) is often too small. A convenient way 
to solve this problem is to introduce a space of test functionals in (L 2 ) and to use its larger 
dual space. 

We like to work with the space of test functions (S) . So we review the standard con- 
struction of (<S) due to [11]. For a more detailed discussion see [5], [10]. Take one system 
of Hilbertian norms |H p | topologizing S (R) which grows sufficiently fast. Then S (R) is 
realized as a projective limit of Hilbert spaces S p (R) : 

«S(R) = f]S p (R) , 

p>0 

where S p (R) denotes the completion of S (R) with respect to | -\ p . Then the space of tempered 
distributions is 

S'(R)= U 5 -p( r ) . 

p>0 

where the dual norm | ■ | topologizes the Hilbert space S- p (R) . 
One convenient choice is 

\f\ p '-=\A p f\ , /€5(R), (1) 

where 

Af(t) = -f'{t) + (t 2 + l) f(t) 

is the Hamiltonian of the harmonic oscillator. Since (L 2 ) is Segal isomorphic to the complex 
symmetric Fock space T(L 2 ) of L 2 (R), we can identify the Fock space T(S P ) with a subspace 
(S) p of (L 2 ) and define the nuclear space 

(5) = n w P . 

Thus we arrive at the Gel'fand triple: 

(s) c (l 2 ) c (sy . 

Elements of the space (S)* are called Hida distributions (or generalized Brownian function- 
als). It is possible to characterize the spaces (S) and (S)* by their S- or T-transforms 
($6(5)*, f€S(R)): 

T<D (/) = «<D, exp (i (-, /)))) = / exp (i (u, /)) $ (cu) dfi (cj) , (2) 

JS'(R) 

S<Z> (/) = (($,: exp (-,/) :)), 

here ((•,•)) denotes the bilinear pairing between (S) and (S)* and we have used the traditional 
notation: 
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:exp (-,/): = C (/) exp ((•,/)),/ € S(R) . (3) 

We denote by E ($) = (($, 1)) the expectation of a Hida distribution <!>. S- and T-transform 
have extensions to the complex Schwartz space 5c (R) and are related by the following 
formula: 

S$(f) = C(f) T$(-if), /G«Sc(R) (4) 

Let us now quote the above mentioned characterization theorem, which is due to Potthoff 
and Streit [16] and has been generalized in various ways (see e.g. [9], [14], [20]). For a full 
proof of a generalized version see [10]. 

Theorem 2.1 

The following statements are equivalent: 

1. F : «S(R) -> C is 

(A) Ray-entire, i.e. for all g, f G S (R) the mapping C3zh F{zf + g) is entire. 

(B) and uniformly of order two, i.e. there exist constants if i, if 2 > such that 

\F{zf)\ <K ie xp(K 2 \z\ 2 \f\ 2 ), feS(R). 

for some continuous norm |-| on S (R) . 

2. F is the S- transform of a unique Hida distribution $ G (5)* . 

3. F is the T- transform of a unique Hida distribution <J> £ (5)* . 

A functional satisfying 1. is usually called a U- functional. 

As an example of an application of this theorem we consider Donsker's delta function. 
Consider the composition 5 a o B (t) of the Dirac distribution 5 a at a £ R with Brownian 
motion B(t), t > 0: 

$ = 5 (B (t) - a) 
$ = <$((-,l [0)t ))-a), a€R. (5) 
The S'-transform of <E> is calculated to be [5]: 



S$ (/) = 




and theorem 2.1 gives immediately that $ is a well defined element in (S)*. 

Now we want to mention some important consequences of theorem 2.1. The first one 
concerns the convergence of sequences of Hida distributions and can be found in [5], [16], 
[10]. 

Theorem 2.2 

Let {F n } n£N denote a sequence of U -functionals with the following properties: 
1. For all f G S (R) , {F n (/)} ng N is a Cauchy sequence, 
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2. There exist K± t K 2 > such that the bound 

\F n (zf)\ < K x exp (K 2 \z\ 2 I/) 2 ) , / G S (R) 

/ioZds for almost all n £ N in a continuous norm \ -\ on S (R). 

Then there is a unique <D G (5)* smc/i £/ia£ T~ 1 F n converges strongly to $. 
This theorem is also valid for 5-transforms. 

Another corollary of theorem 2.1 deals with the integration of Hida distributions which 
depend on an additional parameter (see [5], [8], [10]). 

Theorem 2.3 

Let (Q,B,m) denote a measure space and A 1— > (A) a mapping from to (5)*. Let F (A) 
denote the T-transform of <I> (A) which satisfies the following conditions: 

1. A 1— > F (A, /) is a measurable function for all f £ S (R) , 

2. There exists a continuous norm \ -\ on S (R) suc/i that 

\F (A, < K x (A) exp (tf 2 (A) \z\ 2 |/| 2 ) , / G S (R) 
witfi G L 1 (tt, m) and K 2 G L°° ($7, m) . 
Then <I> is Bochner integrable in some {S)_ q and thus 

$(A)dm(A) G (5)*. 
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T-transform and integration commute 

r(|$(A) dm (A)) (/) = I T($(A))(/)dm(A). 

Again the same theorem holds for the S'-transform. 
Example 

5(B(t) - a) = — [ e iX(B ^ a UX 
2tt Jr 

in the sense of Bochner integration (see, e.g., [5]). 
Remark 

For later use we have to define pointwise products of a Hida distribution $ with a Donsker- 
Delta function 

S ((u,g) - a) 

If the mapping A 1 — ► T$(f + Xg) is integrable on R the following formula may be used to 
define the product $ • 5 

T (cD • 5((u, g) - a)) (/) = ^ J r e~ iXa T«D (/ + Xg) dX, (6) 

in case the right hand integral is indeed a U-functional. 



3 The Feynman Integrand as a Hida Distribution 



We follow [3] and [6] in viewing the Feynman integral as a weighted average over Brownian 
paths. These paths are modeled within the White Noise framework according to 



x (t) = x + J — [ uj (r) dr, 
V m Jt 



in the sequel we set % = m = 1. 

In [3] the (distribution-valued) weight for the free quantum mechanical propagation from 
x (to) = %o to x(t) = x is constructed from a kinetic energy factor exp ^| J t * Q u 2 (r) drj 

and a Donsker delta function S(x (t) — x). Furthermore a factor exp ( ^ J" o a; 2 (r) dr ] is 
introduced to compensate the Gaussian fall-off of the White Noise measure in order to mimic 
Feynman's non-existing "flat" measure D°°x. Thus in [3] the Feynman integrand for the free 
motion reads (the N indicates appropriate normalization): 



J = Nexp ^lil jf* J- ( T ) ^ 6 ( x ^ _ x j _ 



1 

/to 

As in [3] Jq is a Hida distribution, with T- transform given by 



where A = [to,t] and /a, /a c denote the restrictions of / to A and its complement A c 
respectively. Furthermore the Feynman integral E (Jo) = TJo (0) is indeed the free parti- 
cle propagator ^ = = j exp [ 2 | t-t | ( x ~ x o) 2 • Not only the expectation but also the T- 
transform has a physical meaning. By a formal integration by parts 

TJ (/) = E(j e -^M'M*-) e «/(*)-«(>/(*>) e -±I/a< 

The term e Jt o would thus arise from a time-dependent potential W (x, t) = f(r)x. 

And indeed it is straightforward to verify that 

6(t-t )TI (f) = K { f) (x,t\x ,t ) e -/(*)--o/«o) e -i|/Ac| 2 5 (7) 

where 

r M) , ■< 9(t-to) 

exp ^ix f (t ) - ixf (t) - | |/a| 2 + ^ f(r)dr + x- xo^j 

(f) 

is the Green's function corresponding to the potential W, i.e. JT obeys the Schrodinger 



equation 



(idt + ~ f (t)xj K o^ (x,t\x ,t ) = i5(t-t ) 5(x- x ) 
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More generally one calculates 

9(t - to)T (l "ft 1 S(x (U) - x t )) (/) = e-^e te '«- te «"<*°> n fi K W (x^x^t^) . 
V i=i / i=i 

(8) 

Here to < h < ••• < t n < t n+ \ = t and x n+ \ = x . 

In order to pass from the free motion to more general situations, one has to give a rigorous 
definition of the heuristic expression 

I = Iq exp 

In [8] Khandekar and Streit accomplished this by perturbative methods in case V is a 
finite signed Borel measure with compact support. This construction was generalized in [13] 
to a wider class of potentials by allowing time-dependent potentials and a Gaussian fall-off 
instead of a bounded support. 

The starting point is a power series expansion of exp ^— i J t * V (x (r) , r) dr^j using 
V (x (r) , r) = j dxV (x, r) 5 (x (r) - x) : 

exp (-i f V (x (r) , r) dr) = V {-if [ d n t ft / dx t V ( Xi , U) 5 (x (U) - Xi ) 
V Jt / n=0 JA n i=l J 

where A n = {(ti, ...,<„) 1 1 < h < ... < t n < t}. 

In order to consider singular potentials V is no longer taken to be a function V but a 
measure za Under suitable conditions on v it is proven in [8] and [13] that 

/v = / + EH) n / / f n^rftj^on^ife)-*;) 

n=l ^R"^A„ Vj=l / j=l 

exists as a well-defined element of (S)* using theorems 2.2 and 2.3. 



4 The unperturbed harmonic oscillator 

In this section we first review some results of [3] which are necessary for the formulation and 
proof of our main result. Then we prepare a proposition on which we base our perturbative 
method. 

To define the Feynman integrand 

I h = Jo exp (-iJ*U(x(r)) dr^j , U(x) = ^k 2 x 2 

of the harmonic oscillator, at least two things have to be done. 

First we have to justify the pointwise multiplication of Iq with the interaction term and 
secondly it has to be shown that E(ift) solves the Schrodinger equation for the harmonic 
oscillator. Both has been done in [3]. There the T-transform of Ih has been calculated and 
shown to be a [/-Functional. Thus Ih G (<£)*• Later we will use the following modified version 
of their result: 



2sin/c|A| 



Xq + x 2 
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• COS 



k I A| - 2x x + 2x I dt'f (t') cos k (t' - to) - 2x f dt' f (t') cos k (t — t') + 

J to J to 

ft rsi 

+2 / dsi / ds2f (s\) f (S2) cos k (t — s±) cos k (s% — to) 

J t J t 



(9) 

Ho Jto J ) 

7T 

with < k |A| < — , which is easily seen to be a [/-functional. 
For our purposes it is convenient to introduce 

K h (x,t\xo,to) = 0(t-t ) TI h (f) ■ expi |/ A c| 2 • exp (ix f(t ) - ixf(t)) , 

which is the propagator of a particle in a time dependent potential ^k 2 x 2 + xf(t). This allows 
for an independent check on the correctness of the above result. In advanced textbooks of 
quantum mechanics such as [7] the propagator for an harmonic oscillator coupled to a source 
j (forced harmonic oscillator) is worked out. Upon setting j = f their result is easily seen to 
coincide with the formula given above. 

Proceeding exactly as in the free case (see [8], [13]) we first have to define the (pointwise) 
product 

n 

h II 8{B{tj)-Xj) 

3=1 

in (S)* . The expectation of this object can be interpreted as the propagator of a particle in 
a harmonic potential, where the paths all are "pinned" such that B(tj) = Xj , 1 < j < n. 
Following the ideas of the remark at the end of the section 2 we will have to apply (6) 
repeatedly. But due to the form of Tlh(f), which contains / only in the exponent up to 
second order, all these integrals are expected to be Gaussian. 
Using this we arrive at the following 

Proposition 4.1 For xq < Xj < x, 1 < j < n, to < tj < tj+\ < t, 1 < j < n — 1, 
IhYYj=\ ^ (B(tj) — Xj) is a Hida distribution and its T -transform is given by 



T(i h ft s(B( tj )- Xj ))(f) = e -^i/A C iy(^w(t )) n ^ K (f) (x t | X t 
v .7=1 7 3=1 



Proof. For n = 1 we may check the assertion by direct computation using formula (6). 
To perform induction one needs the following 

Lemma 4.2 Let [to,t] C [t ,t] then 

((^ +A V «'])') (f) 

°" W (x , to\x, t) = Kl n (x , t \x, t) , VA G R . 

The Lemma is also proven by a lengthy but straightforward computation. On a formal 
level the assertion of the lemma is obvious as both sides of the equation are solutions of the 
same Schrodinger equation if [to, t] C [t , t ] . □ 

The proposition states what one intuitively expects, ordinary propagation from one in- 
termediate position to the next. 
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5 The Feynman integrand for the harmonic oscillator 
in an external potential 

In this section we construct the Feynman integrand for the harmonic oscillator in an external 
potential V(x,t). Thus we have to define 



ly = ift • exp i^-i V(x(t), r)d 



As for the free particle we introduce the perturbation V via the series expansion of the 
exponential. Hence we have to find conditions for V such that the following object exists in 

(sy 

oo „ „ n 

Iv = h + Y, Hf / d n x / d n t J] V (xj,tj) 5 (x (tj) - Xj ) I h . 

n=l Jnn jA " j=l 

Since we want to study singular time-dependent potentials, we consider v a finite signed Borel 
measure on R x A. Let v x denote the marginal measure 

v x (A G B (R)) = v(A x A) 

and similarly 

u t (BeB (A)) ei/(RxB). 

The following theorem contains conditions under which the Feynman integrand ly exists as 
a Hida distribution. 

Theorem 5.1 Let v = V+ — V- be a finite signed Borel measure on Rx A where the marginal 
measures \v\ x := (u + + v-) x and \v\t satisfy 

i) 3R > 0, Vr > R : \v\ x ({x : \x\ > r}) < exp (-fir 2 ) for some f3 > 0; 

ii) \u\ t has a L°° density. 

Then 

Iy = I h + jr(-i) n [ I ( fl^dx j ,dt j ))l h f[S(x(tj)-x j ) (10) 
is a Hida distribution. 

Remark: Conditions i) and ii) allow for some rather singular potentials, e.g. e ~ n $ n ■ ^ or 
a cut-off interaction, i.e. compactly supported v x , condition i) is of course valid. Note also 
that v is not supposed to be a product measure, hence the time dependence can be more 
intricate than simple multiplication by a function of time. 

Proof. 

1. part: In the first part of the proof we have to perform some technicalities which are 
necessary to establish the central estimate (11). We have to use a very careful procedure to 
achieve that (11) survives n-fold integration and summation in the second part of the proof. 
From proposition 4.1 and the explicit formula 9 we find 

T(h U 6 (B( tj ) - Xj) ) (zf) | < e^'lo (Jn yj^) exp ((|z n+ i| + |x |)| \z\ sup |/|) • 
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cxp I < ^2 ikzxj 



3=1 



sin A: \Aj\ JAj 



dtf (t) cos k(t-t 



j-D 



sink Jaj+i 

n 1 2 

-v 7T \Z\ 



(it/ (i) cos k (t — tj + i) 



We define 



and 



With these 



X = sup |xj| 
0<j<n+l 



= sup|/|+sup|/'| + |/| 

A A 



t(i h ft 5 

V .7=1 7 



exp I <^ ^ ifczxj 



< e 2 



j£l llfll 2 



fn+1 



exp lx*\z\ 11/11 + ™-|A| ||/f 



j'=i 



sin k | A j | j a 



dtf (t) cos k(t-t 



- sin fc |A J+ i| /A J+ ^ /(t)COSMt " tj+l) ]}, 
To estimate the last factor we proceed as follows: 

f dtf (t) cos fc (t - tj) - — ^- / dt/ (t) cos fc (t - t j+1 ) 

|A,-+i| JA j+1 smk |A 1+ i| Ja,+i 



3=1 



sin A; | Aj+i| 



1 /" 

< VfclzlX / dtf(t) / k sin k(t- t) dr 

~ ^ sinA;|A J+1 | 7 A , +1 W it, 



< £>|z|Xsup|/|||A j+1 | 

< |N^||/|||A| 



To obtain a bound for the remaining term 

— — , / dtf (t) cos k(t — tj-i) —3- r / dtf (t) cos k(t — tj 

mklAjlJA, W V J ^ sinfc|A 7 - + i| i A ,+i W V J 



3=1 



we expand = $ dt f(t) cos A; (i — ij-i) and G(tj + i) = f^ +1 dt f(t) cosk (t — tj 

around t,-. This yields with ^ 6 A^ and rjj + \ 6 A J+ i 



< jfe|z|X 



j=i 



| A, 



|A j+ i| 



sinfe|Aj| sin A; | Aj_|_i| 



+ 
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+k\z\Xj2 

3=1 



' ' (/' {Vj+i) cos k {Vj+i ~ tj) - kf (r} j+1 ) sin k (r) j+1 - tj)) 



-k\z\Xj2 

3=1 

Since sup (£ 



For the second term we obtain the bound 

\z\ X 1 -^ sup |/'| + sup |/| + £ | A| sup \f'\ + ^ | A| sup |/| < 

4 a 4 a 4 a 



2sinfc |Aj + i| 
)' = 1 then the first term above is bounded by 
2A;|A| \z\X 



4 a 



< 



\z\X\A\tt 



(ll/ll (2 + fc 2 |A|+A:)) 



Putting all of this together we finally arrive at 
T(l h ft 5 (B(tj) - Xj ) ) (zf) | < exp (L \z\ X ||/|| + g |A| + ±) 



where L = tt + fvrA; |A| + 2k |A| + § | A| (2 + k 2 | A|) is a constant. 
Hence we have the following estimate 



t(i h _ft <^ - xj) ) (z/) | < y nj ) exp ( x * 7 ) exp 



1 7T 



+ -|A| + 



2 2 



2 7y / 

(11)' 
where 7 > 0. 

2. part: In this final step we use the method developed in [13] to control the convergence 
of (10). Although the slight modification to our case is easy we give the basic steps for the 
convenience of the reader. 

In order to apply theorem 2.3 to perform the integration we need to show that 



'n+1 
k3=1 



4|A 3 | 



exp 



{xh) 



is integrable with respect to v. To this end we choose q > 2 and < 7 < ^. With this choice 
of 7 the property i) of v yields that exp (7X 2 ) G L q (R n x A n , \ u\) and with 

1 

Q = J \v\ (dx, dt) exp [^qx 2 ^j 



we have 



/ f ft M {d Xj ,dtj) exp ( 7 ?A 2 ) ) < exp (7 (x 2 + x 2 )) Q n < 

L-/R 11 •'An 7 = 1 / 



OO. 
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Now we choose p such that | + | = 1. Using the property ii) of v and the formula 



n+l 

n 

An j = l 



1 



4\tj-tj-l\ 

we obtain the following bound 



r(i-a)\ n+1 \ A \^-*)-* 



r((n + !)(!-«)) 



a < 1 



n+l 



n ih (dxjidtj) n 



4 |t 7 - - 



< 



< \u\ n/p 



p( 2-p\— | A |2_l( n+ l) 



4^r((n+l)^)^ 



< oo 



\v\ too is shorthand notation for the essential supremum of the L°°-density of \v\t which exists 
due to condition ii). 

Finally an application of Holder's inequality gives 



'n+l 



cxp 



< 



< exp (jx 2 + 7X 2 ) Q n \v 



F (2-p 

n/p V 2 



|A|?^ (n+1) 



too 



2 n+i r ((n + l) 



^— = C n < cc 



Hence theorem 2.3 yields 



R™ JA n j = l 



n * (dxj, dtj) [hUHB (tj) - Xj ) e (sy . 



3=1 



As the C n are rapidly decreasing in n the hypotheses of theorem 2.2 are fulfilled and 
hence 



Iv = E !n e (<S)* 

n=0 



□ 
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